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In many experimental situations, a physical system undergoes stochastic evolution which may 
be described via random maps between two compact spaces. In the current work, we study the 
applicability of large deviations theory to time-averaged quantities which describe such stochastic 
maps, in particular time-averaged currents and density functionals. We derive the large deviations 
principle for these quantities, as well as for global topological currents, and formulate variational, 
thermodynamic relations to establish large deviation properties of the topological currents. We 
illustrate the theory with a nontrivial example of a Heisenberg spin-chain with a topological driving 
of the Wess-Zumino type. The Cramer functional of the topological current is found explicitly in 
the instanton gas regime for the spin-chain model in the weak- noise limit. In the context of the 
Morse theory, we discuss a general reduction of continuous stochastic models with weak noise to 
effective Markov chains describing transitions between stable fixed points. 



I. INTRODUCTION 

Dynamics of complex systems, in the presence of disorder and under external forces, is often modeled by non- 
equilibrium stochastic processes. Despite the considerable interest in such models, exact results, or even effective 
approximation methods are not readily available for generic situations. Instead, specific results were derived under 
restricting assumptions; for instance, the case of steady-states was investigated in many publications (a complete 
list would be prohibitively long), drawin g u pon the methods provided by the Large Deviations Principle (LDP), 
developed, for example, in Refs. H, [2(| I22I |23|. Ii5ll . 

More recently, the LDP-based methods were applied to the case where the quantities of interest are time-averaged 
observables, like currents or densities of particles [g, 0, H, B El, El, El, EL S3 ■ Even more attention was paid to the 
production of entropy, which is actually a linear functional of currents, and related fluctuation theorems. (See, for 
example, Refs. El, EL El El El, El, EL El S3, SI, El, EJ)- However, prior studies have only considered processes 
in spaces with a trivial topological structure, in the sense that we will explain in the following. To the best of our 
knowledge, the topological nature of global currents (fluxes) has not been discussed before. 

In most cases, stochasticity of dynamics and observation errors make detailed knowledge of system trajectories 
unnecessary and distracting. Fortunately, one can often find topological characteristics of the system, which are easier 
to observe. We suggest that properly defined stationary currents are such topological characteristics. 

In the present paper, extending on our recent preprint (l6| . we address the problem of deriving statistical properties 
of empirical (time-averaged observable) currents for non-equilibrium stochastic processes which are equivalent to 
random maps between compact spaces with nontrivial topology. We employ a method inspired by the LDP, but also 
special field-theoretical tools developed originally in the context of nonlinear sigma- models [l[ . 

Our main results concern deriving via LDP explicit asymptotic expressions for joint distribution function of current 
density and density and distribution function of total topological current in non-equilibrium steady state stochastic 
systems (stochastic maps) defined on compact spaces. Existence of a nontrivial large deviations distribution for net 
currents and the resulting thermodynamic relations are intimately related to nontrivial configurations of maps from the 
base domain space to the target space. We mainly focus on yet unexplored relationship between the large-deviations 
probabilistic techniques for non-equilibrium systems and the topological structure of the configuration spaces of the 
model. From that perspective, this article makes a completely novel contribution. In the present study we consider 
a global static violation of detailed balance, in contrast to stochastic pumping problems, where changing parameters 
in time brings the system out of equilibrium 0, 

The power of the aforementioned general results and technique is illustrated on the enabling example of a circular 
spin chain, corresponding to a stochastic process mapping from the torus to the sphere. (See Fig. Q]for illustration.) 
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FIG. 1: Top: consecutive configurations of the string as it wraps around the sphere. Bottom: the same configurations in the 
discretized representation of the string by a cyclic spin chain. 

A physical realization of this model could be the following: consider a nano-structure represented by a circular spin- 
chain of N S> 1 classical interacting spins characterized by unit vectors rii where i = 1, . . . , N, coupled to a stochastic 
external magnetic field B. Note that such a device could be used as a magnetic field detector, by measuring the 
response of the spin chain as a function of time. It is most useful to consider the low-energy, long-wavelength limit of 
the problem, which would correspond to a high sensitivity of the device. Therefore, we will restrict our attention only 
to the spin couplings relevant to the long-wavelength approximation. The spin system can be driven in other ways 
than by an external magnetic field. For example, a nonconservative force, a so-called spin transfer torque, is caused 
by a spin-polarized current [Hi], HH ■ 

Other possible realizations of the spin-chain model are molecular motors [1, [H, |46[ which often exhibit periodic 
motions resulting from nonconservative driving. Although specific examples may not have been yet discovered, we 
believe that there are relatively simple non-equilibrium bio-molecular systems whose functioning is controlled by 
topological currents, probably far more complex than those discussed here in the context of the spin-chain model. 

Performing precise measurements in such setups, in the presence of fluctuations, is essentially related to the ability 
to detect collective modes of the chain rij, by integrating the response over a time interval. Therefore, such a device 
would be a natural detector of the total current associated with the integrated response of the entire chain. Our 
analysis will focus on the probabilistic description of the empirical current (time-integrated response) and generalize 
to compact spaces the study of the current density distribution (or a ID current), the focus of earlier works @, 0, 
H, H, 0, HH, EL EH, 52 m both single- and many-particle systems. The total current in our spin-chain example 
originates from global, topological characteristics of the system, as locally the stationary force driving the system is 
potential. In fact, this example represents the case of the topological driving associated with a multi- valued potential 
of the Wess-Zumino type [l|, H3, [6l|, and thus the essential part of our analysis will be devoted to establishing large 
deviations characteristics of the global (topological) current in the spin-chain model. 

The material in the paper is organized as follows. In Section[TT]we introduce a general stochastic model and a specific 
example with a nontrivial topology. All the results of the paper are briefly discussed and listed at the end of this 
introductory Section. In Section [IIII we present a topological picture of stochastic currents. In Section HVl we consider 
overdamped continuous stochastic processes. We review some results related to statistics of the empirical density, and 
current density, as well as other extensive observables derived from them. We also extend the large deviation theory 
to nontrivial topologies and develop a general approach to obtain large deviation functions for topological currents. 
Section [V] is devoted to a non-equilibrium stochastic model with a nontrivial topology: a classical Hcisenbcrg spin 
chain driven by a topological term of the Wess-Zumino type. In Section IVI1 in addition to giving a summary and 
conclusions, we discuss a general reduction of continuous stochastic models in the weak-noise limit to effective Markov 
chains with the help of the Morse theory. Appendices contain some technical details and auxiliary material. 

II. MODELS AND STATEMENTS OF THE RESULTS 

A. General Stochastic Model 

We consider a stochastic process/trajectory r/ t = {?7(r)|0 < r < t}, or more rigorously 774 : [0, £] — M of duration 
t, which occurs in a configuration space M, i.e., r/(r) € M for < r < t. The configuration space M is assumed to 
be a manifold of dimension m = dimM, so that the particle position t](t) at any given time can be characterized by 
a set rf{r) of local coordinates with i = 1, . . . , m. Our stochastic process can be described by the following Langevin 
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equations 

V i (T)=F*{T,)+?(r l ,T), (1) 

which is a continuous limit of the well defined discrete-time stochastic differential equations written, for example, in 
the Ito form. The quantities discussed below that depend on the continuous time should be understood as limits of 
the their properly discretized forms. In Eq. |T]) F l denotes the deterministic (advection) component of the particle 
velocity, linearly related to the driving force Fj (overdamped dynamics), 

F = //•'/••, , (2) 

via the mobility tensor g %i '(rj) that can be viewed as a Riemann metric on the configuration space M. Due to the 
Einstein relation (fluctuation-dissipation theorem), the same tensor Kg* 3 ', weighted with a factor k that controls the 
noise strength, characterizes the correlations of the Gaussian Markovian noise in Eq. ([3]) : 

(C(r,,r 2 )e (V,n)) = *9 ij (v)S(T2 -ti) . (3) 

If our stochastic dynamics is interpreted as a result of elimination of fast components in harmonic bath modeling 
(to achieve the Markov limit), the Einstein relation means that the bath is at equilibrium at temperature At, and 
non-equilibrium features of the system's stationary state can result only from the non-potential nature of the driving 
force F. Hereafter we imply summation over the repeating indices and assume, without loss of generality, that the 
metric is curvature-free. Eqs. ((3]) are consistent with the stochastic (Onsager-Machlup j4^|) action 

S(Vt) = ±-l drg lk (ff - F( V )) (f, k - F k (v)) (4) 
ZK Jo 

(with giji gijg^ k = 8f) defining the probability measure over Tj t , such that the stochastic average of a functional •(rjt) 
of r\ t (e.g., an observable accumulated over time t) is evaluated according to 

= /M ^( T )exp(-^)) ' (5) 

where the denominator is usually called the partition function. In Eq. ([5]) we use standard notations for the path 
integrals over trajectories and assume proper discretization over the time interval [0, i]. Distinction between different 
discretization conventions is irrelevant for the Cramer functions in the low- noise limit. 

We introduce the empirical density and current at a point x of the trajectory's configuration space: 

Pt(r) U x) = t- 1 J*dT5(x-r)(T)), (6) 
Jt(vt,x) = t- 1 f* drf,*S(x - v(r)). (7) 

We are interested in the large-deviation limit of the joint probability distribution function for p t and J t , 

Vt(J,p) = {S(p t -p)S(J t -J)} c . (8) 

Assuming that the observation time t is large and focusing primarily on statistics of pt, Jt defined above one observes 
that a distinction between an open trajectory with r](0) ^ rj(t) and a closed trajectory with tj(0) = rj(t) disappears 
at t — > oo. This fact, also discussed in detail in Section HTT1 allows us to focus on the analysis of closed trajectories. 

The joint distribution function of current density and density defined in Eq. l[5]) is a very useful and rich object 
carrying sufficient amount of dynamical and topological information about the system one would normally be interested 
in. However, from the point of view of experimentally and computationally desirable low-dimensional characterization 
of the stochastic system, the functional Vt(J,p) is still too complicated. One would like to introduce a version of 
Eq. ([5]) of smaller dimensionality, with inessential parameters integrated out. We suggest that a topologically protected 
quantity satisfying these requirements is the equivalence class components of the intersection index between a closed 
trajectory r\ t and a cross-section a of M , which we denote symbolically as 

^ = / Jt (9) 



Here [a] indicates that the object is invariant with respect to continuous transformation from one cross-section to 
another within the same equivalence class. Obviously, currents can be added and multiplied by numbers, where 
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respective operation is executed over the corresponding current densities. Therefore, currents u) t can be viewed as 
vectors that reside in a certain vector space, a current space. The current components ur t , introduced in Eq. |9|). 
are labeled by linearly independent equivalence classes [a] of cross sections, whose number defines the dimension of 
the current space. Further details (including formal definitions) will be given in Section [IIII As in other parts of this 
manuscript, our main focus will be on evaluating the LDP asymptotic for the respective distribution function 

v t w)^ms(j t a] -J^)) . (io) 

\ lot] I t 

Notice that the number of the linearly independent equivalence classes [a] (the dimension of the current space) 
naturally depends on M, being a small number for common topological problems (field theories). For example, for 
the problem considered in Section [V] the current space is one-dimensional. 



B. Circular Spin Chain 

In the large N limit, N — > oo, we can describe the circular spin-chain model using a map n(y, r) = 
(n a (y, r)|a = 1, 2, 3 & ^2 a n a n a = 1), which represents the three-dimensional unit vector parameterized by the angle 
y £ [0, 2tt] at the time r. By imposing periodic boundary conditions in y and r, we also find convenient to think 
about the model as of a (1 + 1) field theory, i.e. as of a stochastic map S 1 x S 1 — » S 2 . 

To illustrate the general topological results we will consider the following version of the model (UJ : 

d T n = F(ri)+£, F(n) = F c {n) + u[n,d y n], F c (n) = v (d 2 n + (d y n ■ d y n)n) , (11) 

(e(yv,n)e(y2; r 2 )) = k5{ Ti - r 2 )8{ yi - m ) (s ab - n a fa., n )n b fa.,T l )) , (12) 

where parameters u and v correspond to strength of driving and overall noise normalization respectively and [•,•] 
is the standard notation for the vector cross-product. (Eq. (fT2j) guarantees the transversality of the noise term to 
the n(y, t) field, and it is also straightforward to verify that all terms on the rhs of Eq. (fl"2"|) are in fact transversal 
to n(y,t). As argued in the Introduction, this model can describe the long- wavelength limit of a nano-scale spin 
device manipulated by magnetic field [H, 58]. In fact, Eqs. (fTTj) represent the most general weak-noise, weak-driving 
stochastic equations which may be built for the map n(y, t) in the long-wavclcngth limit, thus keeping only low-order 
spatial gradients (over y). It is important to emphasize that the local force in the spin-chain model is conservative, 
and the non-equilibrium character of the stochastic process is due to global (topological) effects. These topological 
aspects of the model will be discussed in detail in the two first Subsections of Section [Vj 

The Onsager-Machlup action (JT| that corresponds to the Langevin dynamics described by Eqs. (fTTj) and (fT2"|) has 
a form 

S(n) = f dr f dy (d T n — u[n,d y n] — v (d 2 n + (d y n ■ d y n)n)) 2 (13) 

Setting v = turns this action into the known model, often called a (1 + 1) nonlinear a- model on a sphere with a 
topological term (often referred to as a 0-term) [l|: 

%)U = -^ / dr ( dy((d T n) 2 +u 2 (d y n) 2 -2u(n-[d T n,d y n})). (14) 
Zk Jo Js 1 

The nonlinear a- model of Eq. (|14p . considered as an imaginary-time field theory, has ultraviolet divergences, and 
thus requires a small scale regularization. In fact, our circular spin chain model (|11|) may be viewed as a regularized 
counterpart of the u-model, where the number N of spins plays the role of the ultraviolet cut-off parameter. An 
intuitive explanation for this regularization is that the v-tevm in Eq. (|lip acts to align all the spins and, therefore, 
suppresses the short-range fluctuations. 

In this manuscript we will focus on analysis of the circular spin-chain model in the wcak-noisc limit k <C \u\ < v, 
which will also be coined (by the reason to be spelled out later in Scction|V| the limit of "instanton gas" , correspondent 
to moderate topological driving . We demonstrate that the configuration space of the spin-chain model has one 
topologically nontrivial cycle (i.e., the current to is single-component). In the considered limit the system spends most 
of the time around its stable configuration, whereas the current is generated by rare events, referred to as instantons, 
whose interaction can be neglected due to long time intervals between them. 



1 Another "anomalous" regime of the sigma-model, v <C K <C |tt|, was analyzed by Polyakov and Wicgman [49l. |50H in the zero topological 
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C. Statement of Results 

The main results reported in this manuscript are as follows: 

• In Section UlTl we establish a topological nature of the average current generated over a long time in a stochastic 
system. The topological stochastic current u> resides in a vector space, referred to as the current space, whose 
dimension is given by the number of independent 1-dimcnsional cycles of the system configuration space M. We 
demonstrate that there are two equivalent ways to view the generated topological currents: (i) the rates with 
which the stochastic trajectory loops around the independent 1-cycles, and (ii) the equivalence classes of the 
divergence- free current density distributions. The equivalence of these two views is established by the Poincare 
duality represented by Eq. (fT5|) . 

• We show in Section HVl that for a stationary stochastic process the joint probability distribution for the empirical 
current density and density in the t — > oo limit takes the form Vt(J,p) ~ exp(— tS( J, p)) with the Cramer 
functional 

s(j, P) =[ (15) 

Jm 2k -p 

This generalizes the previously reported results @, 0, H, H EH, El, E3] to the case of topologically nontrivial 
compact spaces. 

• In Section HV CI we describe a general method to derive large-deviation statistics of particular currents (e.g., the 
Cramer function of the topological currents) from Eq. (|15[) in a variational way by solving a respective set of 
equations for currents and densities. 

• We illustrate the utility of the general approach on the example of the thermodynamic limit N — > oo of the 
spin-chain model defined above in Section III Bl We show that the topological current space of the model is 
one-dimensional, i.e., the generated current is described by uj € R. We compute the Cramer function S(u>) in 
the weak-noise limit n <C v, \u\ for not too large values of the generated current uo for \w\ <C |«|(ln(|i;|/«;)) _1 , 
when the instanton gas mechanism dominates the current generation. In this instanton-gas regime the Cramer 
function S(uS) has the same form as for the effective ID random walk (equivalent to a circular two-channel 
single-state Markov chain) with jump rates k + and k_ in opposite directions: 

Cf \ i ,1 ° J+ V / ^ 2 + 4k + K - / 2 i A MR\ 

o \UJ) = k+ + K_ + OJ In y u + 4k + k_ , (lb) 

the rates k± being expressed through the parameters u and v of the model (fTTj) by means of Eqs. ([62]) . (1891) . 
(EU), HMD, and 

• In general, the vector of topological currents u) is not related to the work produced by the driving force, 
f Q d,T(F(r)) ■ rj). However, in the case of topological driving, as in the spin-chain model, the work becomes a 
linear functional of the current vector. 

III. TOPOLOGICAL VIEW OF STOCHASTIC CURRENTS 

Historically, a concept of current appeared in physics on a macroscopic level as a way to describe a flux of any 
kind. Consider an electric circuit, represented by a circular wire with a static electric field (for example, provided by 



charge sector (not relevant for our application) . The authors mapped the model onto a system of infinite-component massless interacting 
fcrmions, analyzed it with the Bethe ansatz approach and thus arrived at a remarkable exact and nontrivial solution. Notice, that an 
intermediate case of "interacting instantons", correspondent to strong topological driving and vanishing noise, k < « < |ti|, constitutes 
yet another interesting regime, which to the best of our knowledge was not studied yet. The Cramer function in this later case can 
be calculated explicitly by integrating over the instanton solutions and Gaussian fluctuations around them in the limit v = 0. As 
demonstrated in Ref. 12811 in the context of the cr-model as a field theory, the problem is equivalent to finding the ground-state energy of 
the corresponding (1 + 1) sine-Gordon model, where v will play a role of the ultraviolet cut-off parameter. In the case v < u the saddle 
and stable points of the spin-chain system coalesce producing a peculiar fixed point at 8 = which is stable when approached from one 
direction and unstable from the other one. All such situations are beyond the scope of this work. 
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a battery) where the electrons on average move in one direction. The current is defined as a charge crossing some 
oriented section a per unit time. The current can be considered as a sum of currents from individual particles. An 
individual contribution uj is given by ui ~ Nt" 1 with TV = N + — 7V_, where N± is the number of times the particle 
trajectory rj crosses the section a in the positive and negative directions, respectively (see Fig. 




FIG. 2: A single open trajectory in the circular wire. 

The configuration space M for the particles in the circular wire is three-dimensional, i.e., m = dimM = 3, and can 
be represented as M = S 1 x D 2 (with S 1 and D 2 being a circle and a two-dimensional disc, respectively), whose 
boundary is dM = S 1 x S 1 . If the particle trajectory rj is closed, i.e., r] : S 1 — » M, the number TV does not change 
upon deformations of the trajectory and the cross-section a. This number, known as the intersection index, only 
depends on the equivalence classes [rj] and [a] and can be denoted by [77] * [a] 6 Z. The equivalence, based on 
deformations, is the homotopical equivalence 2 . 

A. Intersection Index and Stochastic Currents 

The intersection index does not change if any of the cycles is replaced by a homologically equivalent counterpart. 
Since cycles can be added and multiplied by the integers (by forming disjoint unions and changing orientations), 
this equivalence can be understood if we define a zero cycle. Formally, a j-dimcnsional cycle is called homologically 
equivalent to zero, if it is a boundary of a (j + l)-dimcnsional region mapped into M. Obviously homotopy equivalent 
cycles are homologically equivalent. The set (actually an Abelian group) of homological classes of j-dimcnsional 
cycles in M is called the j-th homology of M and denoted by Hi(M; Z). In the circular wire, an element of H\{M\ Z) 
represents the number of times the particle trajectory moves around the circuit: i/i(A/;Z) = Z. The current u> 
associated with the trajectory r\ can be naturally defined as u> = [ry]< _1 € Hi(M; R) = Hi(M; Z) (g>zR. In the circular 
wire of Fig. [2] the current has only one component because H±(M;M.) = R. 

The topological picture, presented above on the simple example of a circular circuit, can be extended to a much 
less intuitive general case in a pretty straightforward way by viewing an averaged current uj generated in stochastic 
dynamics in the configuration space M as an element in H\(M; R), the first homology group of M with real coefficients, 
associated with the homology class [77] of a stochastic trajectory. Since H\{M] R) is a real and in most relevant cases 
finite-dimensional vector space, the generated current can be viewed as a vector. The homology groups can be 
finite-dimensional and computable even in the field theory when the configuration space M is represented by an 
infinite-dimensional (in the non-regularized continuous limit) space of maps. In particular, in our enabling case of 
71/ = Map(5 1 , S 2 ) discussed in SectionfV] the current also has only one component since 7Ji(Map(5' 1 , S 2 ); R) = R. 

Notice also that, originating from counting, stochastic currents can be viewed as topologically protected observ- 
ables, available in single-molecule measurements, which provide stable and at the same time rich information on the 
underlying stochastic processes. 

The topological picture of stochastic currents can be formulated in the simplest way when defined for closed 
trajectories (loops). However, a general stochastic trajectory is open, i.e., its end point is typically different from the 
starting one and thus adopting the topological language of intersection indexes and currents may seem problematic. 
In the long-time limit considered in this manuscript, extension from open to close trajectories does not constitute a 
problem. Indeed, when counting the intersection index TV of a long trajectory, one can always close it with a segment 



2 Here and below [■] is used as a notation for the equivalence class of 
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(e.g., a geodesic line) that is much shorter than the trajectory itself, which creates an uncertainty no more than one 
in a big number N 3> 1. In contrast to a boundary effects with unbounded fluctuations [52], [6(J, in the long-time limit 
in our case one can safely ignore changes in the statistical properties of the trajectories caused by this modification. 

B. Currents, Current Densities, Zero-Curvature Vector Potentials and Poincare duality 

Currents can also be represented by the current densities J defined as time integrals over trajectories as in Eq. |7|). 
In this Subsection we will discuss some important "static" relations between J, the current oj, and the intersection 
invariant [a]. 

The time-averaged current density from Eq. ([7]) is a random variable on the space of stochastic trajectories tj(t) 
and satisfies the relation 

div x J t ()7t,x) = r 1 (S(x - vt(t)) - 6(x - r?t(0))) . (17) 

Therefore, our first observation is that the random variable divJ t ~ t~~ x vanishes in the limit t — > oo. Moreover, 
Eq. (fTT|) guarantees that the current density is exactly divergence-free, div x Jt(T)t, x ) = 0, if only closed stochastic 
trajectories are considered, which, as argued in Section IIII Al does not affect the long-time behavior of the relevant 
distributions. 

The vector field of the current density on the m-dimcnsional phase space can be naturally viewed as a differential 
form of rank (m— 1) because its integration over a (m— l)-dimensional cross-section section a results in a current. In 
the following we will use the same notations for vector fields and the corresponding differential forms. The divergence 
of the current density J is represented by the exterior derivative as div x J = dJ, and Eq. (|17p suggests that the 
current density is a closed from: dJ = 0. Therefore, it represents a class [J] in the de Rham cohomology H m ~ 1 (M; R). 
The relation between the homological and current-density representations of the current is determined by the Poincare 
duality i?i(M;R) = H m - 1 (M;R) [5f|. This can be formally expressed as 



V[a] G F m _i(M;Z) : u * [a] = ^[q] * [a 




thus representing that per unit time intersection index t — 1 [77] * [a] of a closed trajectory r) with a cross-section a is 
equal to the integral over a of the current density J produced by the trajectory. 

We conclude this Section by noting that the Poincare duality if m _i(M;R) = H 1 (M;1&) in the complementary 
dimension leads to a natural representation of the sections a in terms of the vector potentials A, which are curvature- 
free, i.e. 3 

(diAj - d 3 Ai)dx l A dx J = 0, (19) 

since the cross-section a represents a homology class [a] G H m -i(M;M), whereas the corresponding vector potential 
A represents a gauge equivalence class [A] G i? 1 (M;K). Here in Eq. (fl9|) we used standard wedge-product notations 
for the differential 1-forms [3^|. The relation "[A] correspond to [a] via the Poincare duality" can be conveniently 
expressed in a way extending Eq. (jTSJ , 

V[J]eH m - 1 (M;Z)\divJ = 0; [j=[ AAJ=[ dxA-J, (20) 

which implies that the current components u> * [a] can be labeled by the gauge classes of the curvature-free vector 
potentials. This expresses the general Poincare duality Hj(M;M.) = H m ~i{M; R) in dimension j = m — 1, i.e., 
if OT _i(M;R) = /^(MjR). Using Eq. it can be formulated as follows, (i) With any (m - l)-cycle a wc can 
associate a closed 1-form A so that Eq. (|2"0|) holds for any divergence-free current density distribution J. (ii) Any 
two forms A and A' that satisfy the condition (i) are equivalent, [A] = [A 1 ], (iii) Homologically equivalent cycles 
generate equivalent forms, i.e., if A and A' are respectively generated by a and a' , then [a] = [a'] implies [A] = [A']. 
Note that the conditions (i)-(iii) define a linear map H m -i(M;M.) — * i? 1 (A/;R). The Poincare duality also means 
that this map is an isomorphism. Note that Eq. (|20|) is noticeably distinct from its counterpart Eq. (| 18|) . The latter 
describes the Poincare duality in the complementary dimension j = 1, by associating the currents ui G H\{M; R) with 
the equivalence classes [J] G _ff m-1 (M;R) of divergence- free current densities J. 



3 Where here and below di is our shortcut notation for d Xi ■ 
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IV. LDP FOR EMPIRICAL CURRENTS 

In this Section, we describe results concerning the Large Deviations Principle (LDP) for (possibly) dependent 
sequences of random variables, including its application to the case of time-averaged current densities, and present 
a path-integral derivation for the relevant functionals. The material of this Section is organized as follows. In 
subsection IIV Al we very briefly review the foundations that stand behind the LDP and formulate a variational 
principle relating the LDP for current density and density to the LDP for conjugated vector and scalar potentials. 
In Subsection ITVBl we describe a convenient representation of the LDP in the intuitive path-integral language. In 
Subsection IIV CI we present a general method to derive the Cramer function of the topological currents from the 
Cramer functional of the joint density and current density distribution. 

A. LDP for empirical currents and the Gartner-Ellis theorem 

Gartner- Ellis (G-E) theorem [26|, H(| formalizes the Large Deviation Principe (LDP). It provides a convenient 
theoretical tool for studying the long-time behavior of stationary driven systems. In this Subsection we discuss 
fundamental relations and objects associated with the G-E theorem and the LDP. 

Consider a sequence of random variables £ T e 7i with r = 1, 2, . . . , t. Denote by Pt(f) the distribution of the average 
ipt = tr 1 Yl T =i ^ t on t ne vector space H. Define the corresponding generating function Qt(ip) = (exp (tip ■ (ft)}, with 
ip E H* being linear functionals in H. If the limit 

AO/0 = I™ i" 1 In (Q f (</>)) (21) 

t — »-oo 

does exist, is represented by a convex and bounded from below function, and 

£(¥>)= sup (V-V-A(VO) (22) 
if- 

is represented by a bounded from above lower semi-continuous function with compact level sets, then {ipt} satisfies 
the Large Deviations Principle (LDP) with the rate (or Cramer) function S. This is known as the Gartner-Ellis (G-E) 
theorem [26|, The LDP is formulated in terms of the probabilities Pt(K) = J K dPt((f) for the random variable 
{ft} to belong to K C H. In many cases it can be formulated in a stronger, yet simpler form (see [H, [3(j for a 
general formulation): for any (measurable) set K C 7i 

lim r 1 In (P t (K)) = S(K), (23) 

t — *oo 

where by definition the Cramer function is 

S(K) = inf S{(p), VK C H. (24) 

Note that even though the LDP and the G-E theorem are formulated above for a discrete set of random variables, 
the formulation also implies straightforward extension to the continuous parameterization t G R + = {t G K|t > 0}, 
where t plays the role of time. In particular, we choose ipt = (pt,Jt) G TL with div«7 t = 0; the empirical density 
and current density are introduced in Eqs. ((6]) and ((TJ), respectively. In accordance with discussion of Section [lll| the 
argument ip = (V, A) of the generating function is represented by a potential function V conjugate to the density, and 
by a gauge equivalence class of an Abelian zero curvature gauge field A conjugate to the current density. The LDP 
can be viewed as a mathematically correct way to formulate a physically intuitive statement that under the described 
conditions at long enough times t the probability distribution adopts an asymptotic form 

P t {ip) - exp(-*5( V )). (25) 

By its formulation, specifically due to Eq. (pM)) . the LDP allows the rate functions for reduced variables to be 
obtained via a variational principle. Let Ti' be the residence vector space for the reduced variables tp' t with the 
reduction (projection) map p : H — > HI . In our case the reduced variables are the currents tp' t = u>t, so that 
TL 1 = H\(M) R) and the reduction map, obviously defined by p(p t ,Jt) = [Jt], is linear. The conjugate to (p' t argument 
if)' of the generating function Q'(ip') resides in (Tt')* = ff 1 (M;IR) and, therefore, can be represented ip' = [A] by a 
gauge equivalence class of a curvature-free dA = (l/2)(diAj — djAi)dx l A dx^ = vector field. Applying the LDP we 
obtain 



S'{u>) = 5(p _1 (w)) = inf S(p, J). 

[J]=UJ 



(26) 



9 



Moreover, the LDP can be further interpreted as effectively implementing Legcndre-type transformations between 
thermodynamic potentials (effective actions). We describe these general relations in the following paragraphs. 

Let us also note that our stochastic theory of maps S 1 x Y — > X is, in fact, a field theory; however, being identified 
as a theory of stochastic trajectories S 1 — > M, with M = Map(Y, X), it can be interpreted as classical one-dimensional 
statistical mechanics in a circular system of the size t with the target space M. The time-averaged current Uf = t Qt 
can be viewed as the density of the topological charge Qt = [SSt(A)/SA], obtained from the variational derivative 
of the gauge-invariant action St(A) with respect to the stationary vector potential A, followed by switching to 
(homology) equivalence classes, the latter operation denoted by square brackets. Naturally, a gauge equivalence class 
[A] G i? 1 (M;R) of curvature free vector field A can be interpreted as the chemical potential that corresponds to 
the topological charge Q. According to Eqs. (|2ip. (|2"5|) the rate S(u>) = i _1 f2 t (o;) and the logarithmic generation 
X(A) = t~ 1 Wt([A]) functions can be interpreted as the densities of the free energy and thermodynamic potential, 
respectively, in the thermodynamic limit t —> oo. Naturally, they are connected via the Legendre transformation 

dn = Sdt + [A] dQ, dW = Sdt-Q- d[A] , (27) 

and the thermodynamic potential is represented by the effective action: 

-W t (A) = / p tS M Jt-A\ = f Vrte -S(r,)+tf M J t -A^ (2g) 



Therefore, applying the LDP argument to the empirical current u t will lead to a rate function S, which is simply 
the generation rate of the effective action for the currents of the theory. This restatement of the problem would be 
rather trivial, unless the theory had interesting topological structure. Indeed, as known from quantum field theories 
in such cases currents may have non-perturbative, anomalous terms arising from global topological effects. As we 
will show in the remainder of the paper, this also happens in the non-equilibrium stochastic theory with driving. 

B. Path-integral picture of LDP and the current density functional 

In this Subsection we, first, derive the general (compact spaces) LDP expression Eq. (fT5|) for the distribution of 
density current and density, and, second, discuss respective transformation to the conjugated variables (potentials). 

1. Derivation of Eq. f 15\) 

Using a standard representation for the Dirac ^-functional in Eq. (jSJ) gives for the probability distribution function 

V t (J,p) ~ JvAVVe- u f dx( - A J+Vp) Jvrj t e- s{r >" A ' V \ (29) 

S( Vt ; A, V) = s( m ) -iJ Q dT (rf t Mm) + v(m)) , (30) 

where the integration is performed over real auxiliary fields A and V, and T>A and W are the standard field-theoretical 
notations for functional differentials/measures. 

In the large deviation limit (t — > oo), the path integral in ([30]) is estimated as 

Vr)te Sim;A,v) = Tre tc A , v _ exp (j;X(A, V)) , (31) 

with — X(A, V) being the lowest eigenvalue of the operator — Ca,v, 

Ca,vR = Ap, where t A y = (n/2)VjVj - VjFj + iV, V, i), (32) 
For the normalized ground state eigenfunction p(x) we obtain 

X(A, V) = J dx£.A,vp(x) and / dxp(x) = I . (33) 

Applying further the saddle-point approximation to the functional integral in Eq. (|29[) with respect to A, V, and 
using Eqs. ([30"|) - (|33"|) . we arrive at the following equations: 

J(x) = (F + iKA-(n/2)d)p(x), p = p. (34) 



10 



Actually, this saddle-point approximation involves a deformation of the integration contours to the complex plane, 
which makes A imaginary and V = in the saddle point. Thus, the saddle-point approximation corresponds to the 
supremum with respect to (— iA) in the G-E theorem (in subsection IIV Al the definition of A differs from that of this 
subsection by a factor i). Solving Eqs. ([551 154"]) for A, V and p and substituting the result back in the saddle-point 
expression for the integral in Eq. ([29]) yields Eq. ([15)) for the Cramer functional. Note that boundary (surface) terms 
do not contribute due to the compactness of the target and base manifolds. In deriving Eq. (|31|) from Eq. (|30|) . the 
gauge freedom was fixed by the requirement that the left cigenfunction of Ca,v, conjugated to the right eigenfunction 
p, equals unity. 



2. LDP for charges generated by scalar and vector potentials 

The explicit large deviation result can be immediately used to get thermodynamics-like relations for Cramer 
functions of derived objects. One introduces the sets A^ a \x) and V^ b \x) of vector and scalar potentials, respectively 
(which can also be interpreted as gauge fields, i.e. generators of continuous symmetry transformations mentioned 
above), and the corresponding sets WA(Jt) and uy(pt) of charges 

w^(J t ) = [ dxA ( f ) (x)Ji(r,;x)=t- 1 [ drrf A[ a) (r,), 
Jm Jo 

v$\p t ) = I dxV^(x) Pt (r 1 -x)=t- 1 f drV^(v(r)). (35) 



M 



At t — > oo, the joint p.d.f. Vt(w, u) = (8(w — WA(Jt))S(u — uv(pi))}$ of WA(Jt) and uv(pt) has the large deviation 
form Va,v(w, u) ~ cxp(— tSA,v(w, u)), where 

Sa,v{w,u)= inf S{J,p). (36) 

wa(J)=w, uv(p)—u 

In the path-integral terms, the variational principe expression (|36[) can be obtained by representing the probability 
distribution 

Va,v(w, u) ~ J VJVpS(w - w A (J))S(u - u v {p))e- ts{J ^ , (37) 

followed by computing the integral in the t — > oo limit using the saddle-point approximation. 

Considering a marginalized version of Eq. (|36p . associated with the distribution functions of the charges, generated 
by the vector potentials only, we have Va(w) ~ exp(— tSA{w)) with 

S A {w)= inf S(J,p). (38) 

W A { J) = uu 



The variational principle, represented by Eq. (|38p has two important implications that correspond to two specific 
choices of the gauge field sets A^ a \x). The choice A(x) = F(x) leads to the observable WF(Jt) that, obviously, 
represents the work (entropy) production rate, whose Cramer function S(wf) satisfies the fluctuation theorem. 

However, in this manuscript we focus mainly on the other implication of Eq. ([55]) associated with the set {A^ a ^} 
of curvature-free dA^ = 0, Va is chosen in a way so that the corresponding set {[A'")]} of equivalence classes forms 
a basis set in i/ 1 (M;R). According to the Poincare duality, as described at the end of section [III Bl and specifically 
due to Eqs. (|T5|) and ([20")) the current lj = [J] as the homology class of the current density J is fully characterized by 
the values of the set w A w (J) of observablcs. Therefore, for the described choice of {A^ a ^} the variational principle 
of Eq. ([55[) is equivalent to the variational principle of Eq. (|26D . 

We should note for completeness that since us £ Hi(M;M), the Fourier variable xp conjugate to u> resides in the 
space ip £ (Hi(M;M.))* = H 1 (M;M.), which results in the following representation: 

P,(w) = (^-w)) r / dipe u ^-^K (39) 



Going along the lines of derivation of Eq. (|30[) and using the Poincare duality ([20|) . which allows the identification 
lj = [J] and i\) — [A] , we recast Eq. ([3"9"|) as 



P t ([J])~ / d[A]e- lt ^M dxA - J f Vrj t e- S ^ A ^ 

JH 1 (M;m J 



H 1 (M:l 



d[A]e- U S*t d*A-J e t\(Afi) ^ 



11 



This representation implies that A ([A]) = X(A, 0) is obtained from X(A, V) by restricting the to the zero V = scalar 
potentials and curvature- free dA = vector potentials. 



C. Derivation of the Cramer functional for topological currents 

This Subsection describes a general strategy for calculating the Cramer functional S{u>) of the topologically pro- 
tected currents u>. This is achieved via the variational procedure, formulated in subsection IIV Al [Eq. ([26]) ]. equivalent 
to "integrating out" current density J and density p dependence for a fixed value of u. 

Specifically, we will minimize the Cramer functional (|15[) S(J,p) over J and p under the following conditions 4 : 

divJ = d t J = 0, Jdxp=l, [J]=u>. (41) 

A variation of the current density that satisfies the continuity condition and keeps the topological current constant 
has a form 8J = d^C with £ being a 2-form. A straightforward calculation allows the requirement SS/6£ = to be 
represented in a form: 

dA = dxA = with A = p~ 1 (Fp- J-{n/2)dp), (42) 

where the second equality in Eq. (|4"2"|) should be viewed as the definition of the vector field (1-form) A. The stationary 
current that minimizes S{u>) corresponds to A = 0. The vector potential A determines how the density and current 
density distributions locally differ at the given topological current and in the stationary regime. Note that the vector 
field introduced in Eq. ([42]) differs from its counterpart introduced earlier in Eq. (j34|) by a factor {—in). Thus, with 
a minimal abuse, we use the same notation for both and hereafter stick to the one given by Eq. (|42| . 

Variation of S{J, p) with respect to p is performed in a straightforward way by introducing a Lagrangian multiplier 
A to satisfy the normalization condition [the second relation in Eq. (|4Tj)]. This results in 

{ti/2)divA + F ■ A {1/2) A 2 = -kX. (43) 

Representing A gradient 

A = —ndlnp-, (44) 

where p~{x) is a newly introduced scalar function, and substituting Eq. (|44[> into Eq. (|43|) . one finds that the quadratic 
first-order Riccati-type equation (|43p is transformed into the following linear second-order differential equation 



C*p-(x) = Xp-{x) (45) 



O = { K /2)8 2 + F-d. (46) 



with the adjoint Fokkcr-Planck operator 



On a compact manifold M, eigenvalues A are discrete and bounded from below. A solution of Eq. (|45|) corresponding 
to the lowest — A determines p~{x) and A. However, the function p~{x) is not necessarily single- valued and can acquire 
uncertainty in the result of going around topologically nontrivial cycles. We can find a unique solution A of Eq. (|46|) 
with the minimal eigenvalue if we fix a set Z = {Zk\k = 1, • • • ,n) of topological parameters defined by the integrals 



\nZ k = n- L / A t {x)dx l (47) 

over the set {sfc}fc=i,...,« of the topologically independent 1-cycles of M. 

To summarize, Eqs. (|4"Tj) -([4"3" |) constitute a complete set of equations that can be used to obtain the Cramer functional 
S(lu) as well as the distributions p and J. This is achieved in three steps. 



4 Throughout the paper we use three equivalent representations of the current density: vector field, 1-form and (m — l)-form. The vector 
field is related to the 1-form through the natural metric tensor that characterizes noise correlations. The relation between the Hodge 
dual differential forms of degrees 1 and (m — 1) is also determined by the metric tensor. 
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(1) We first solve Eq. ((43)) together with the constraint dA = (i.e., the first relation in Eq. (|42[)) with respect to 
the vector potential A. As explained above this is equivalent to solving the linear lowest eigen-value problem 
(|45| with the additional topological freedom fixed unambiguously selecting the sets Z = (Zk\k = 1, • • ■ , n). 



(2) We combine the first two relations in Eq. (jiTj) with the second relation in Eq. (|42]l . which results in 

(n/2)dp — (F — A)p = J, rftj = 0, and [dxp=l. (48) 



^From this linear system we find p and J in terms of Z parameterizing A 5 . 

(3) The Cramer functional S can be obtained as a function of Z by substituting the solution obtained on steps (1) 
and (2) into Eq. (fTS)) . Finally, the substitution of the obtained solution for J into the third relation in Eq. (|4"Tj) 
establishes a relation between Z and a>, which, being resolved with respect to Z in terms of u>, results in S(u>). 

In the case dF = of topological driving, the system of equations Eqs. (|4"Tj) - (|4"3")) for the Cramer functional can be 
further simplified. One applies the operator (n/2)d — (F — A) A to the second relation in Eq. (|4"2")l and makes use of 
the relations d 2 = 0, dA = 0, and dF = 0. This results in ((n/2)d — (F — A)A)J = 0. Combined with the generic 
expression for the Cramer functional (|15|) and applying some reordering of terms, this allows the system of equations 
(gT])-([13]) and Eq. (TJSD to be recast as 

(n/2)divA + F ■ A — (1/2) A 2 = — kA, k' 1 / A l {x)dx l = In Z kl dA = 0, (49) 

J s k 

(n/2)dJ-FA J + AA J = 0, d^J = 0, (50) 

(n/2)dp- Fp + Ap= -J, dxp(x) = l, (51) 

Jm 



S(u) = {2k)- 1 dxpA 2 , / J = w fc , (52) 

where for k = 1, . . . ,n one denotes by Sfc and a k the dual sets of topologically independent 1- and (to — l)-cycles, 
respectively, i.e., [sk] and [ak] form the basis sets of H\(M) and iJ, n _i(il/), with the intersection property * ak> = 
5kk>- Note that in Eqs. (|5TJ|) the current density is understood as the 1-form. 

Therefore, the procedure of finding S(lj) in the case of the topological driving dF = can be summarized as follows: 

(i) For an arbitrary set Z = (Zk\k = 1, . . . , n) we find a unique solution of Eq. (|49p that corresponds to the minimal 
value of A and express the vector potential A in terms of the parameter set Z . 

(ii) Upon substitution of A found on the first step into Eq. ([50]) . the latter can be viewed as a system of homogeneous 
linear equations for the current density J, whose solution is unique up to a multiplicative factor. 

(iii) We substitute the obtained current density J and vector potential A into Eq. (|5ip . which is now viewed as a 
linear equation on the density p with an inhomogcncous term represented by J. Therefore, the normalization 
condition (the second equality in Eq. (|5ip ) determines the prefactor in the current density, which identifies 
normalization (multiplicative) factors both for J and p. 

(iv) We substitute A, J, and p into Eq. ([52"]) . This expresses the Cramer function S and the topological current u> 
in terms of Z . Expressing Z in terms of uj gives the desired Cramer function S(us). 

It is also instructive to note that an alternative representation for the Cramer functions 

S = — A — k ^ j dxJ ■ A = (2k)" 1 J dxpA 2 , (53) 

allows the eigenvalue A to be determined, which is useful for the weak-noise calculation. 

The formal scheme described above will be implemented explicitly on our enabling example of the topologically 
driven spin chain in Section [V CI 



5 This derivation is similar to the one given in Ref. [39l | for the stationary current caused by a force field; in our case the force field is 
replaced by F — A. 
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V. NON-EQUILIBRIUM CYCLIC SPIN-CHAIN 

This section focuses on the applications of the general formalism developed above to the model of topologically 
driven system of N classical spins, arranged in a circular chain. In the thermodynamic. N — > oo, limit the configuration 
space becomes infinite-dimensional, and the system can be viewed as a (1 + 1) stochastic field theory with the target 
space S 2 (sec Fig. [T]). The model has three parameters v, u, and k that describe the relaxation rate, the rate of 
topological driving and the noise strength, respectively. As briefly discussed above, the relaxation term suppresses 
the short-range fluctuations, and, therefore, makes the model divergence-free in the thermodynamic limit. We arc 
considering the weak-noise limit, k <C v, \u\, whereas v and \u\ can be comparable although such that \u\/v is not 
too large. The last requirement translates into the condition that the constant loop n(y) = tiq solution constitutes a 
stable stationary point of the deterministic (zero noise) dynamics. 

This Section is organized as follows. Subsection IV Al is devoted to formulation of the topological driving in terms 
of a multi-valued potential, referred to as the Wess-Zumino potential. We also briefly discuss the finite-dimensional 
(regularized) approximations for the infinite-dimensional field-theory configuration. Some details on the spin-chain 
regularizations are presented in Appendix where we argue that, starting with a large enough N, the relevant 
topological properties of the finite-dimensional approximations stabilize to their continuous limit counterpart. In 
Subsection IV Bl we describe the instanton (optimal fluctuation) mechanism of the current generation and identify the 
structure of the instanton space Mo C M that consists of all the configurations the optimal trajectories (instantons) 
pass through. Subsection IV CI contains the derivation of the main results concerning the spin-chain model. Here, we 
calculate the Cramer function of the topological current u>, generated in the system, by implementing the general 
procedure outlined in Section [IV CI This is achieved by solving the relevant Fokkcr-Planck type equations equations, 
using an ansatz for the density p, current density J, and the vector potential A. The implemented ansatz is asymp- 
totically exact in the weak-noise limit. The Cramer function S(co) is derived for \uj\ <C |u|(ln(|«|/«;)) _ , i.e., in the 
instanton gas regime. 

A. Wess-Zumino interpretation of the Circular Spin Chain Model 

Generally, a system described by Eq. (TTJ is not globally driven if the force field is given by an exact differential, i.e., 
F = —dV with (dV)i = diV, where V is some scalar potential function. Although the force field form is not exact in 
the spin-chain model, it is still closed: dF = 0, with (dF),j = (l/2)(diFj — djFi), or, in other words, the force has 
zero curvature. Therefore, the deterministic force in Eq. (|11[) can be expressed as 

(54) 

n 2 = l 

in terms of a multi- valued potential Vwz{n). Without a single- valued potential, the system is globally driven, and 
its stationary state can only be a non-equilibrium one, with a current being generated. 

To determine the potential, consider some reference configuration n £ M, e.g., a constant loop (a set of collinear 
spins) no(y) = tiq and for an arbitrary configuration n <G M represented by n(y) consider a path x '■ [0^] ~~ * 
represented by x{Vi s ) with \x{y, s)\ = 1 along the path that connects n to n, i.e., x(0) = n and x(t) = n : or > 
cquivalently, x(y, 0) = Uq and x(y, t) = n(y). In the following we will skip the dependence on (y, r) when obvious. 
To derive the potential, we can simply integrate the force F along the path %. According to the Stokes theorem, 
dF = guarantees that we obtain the same potential if the paths x an d x' used in the integration are topologically 
equivalent. Thus, we obtain the potential 

Vwz(n) = V c (n) - uip B (n), V c (n) = - dy(d y n) 2 , (55) 

1 Js 1 

<p B (n) = - dr dy{x ■ [d T x,d y x)] ■ (56) 
Jo Js 1 

The "elastic" globally potential term V c (n), whose variation produces F c (n) in Eq. (JTTJ) , enforces relaxation of the 
spin system to a y-uniform distribution (a constant loop). The second term oc tps is topological and similar to the 
multi- valued Wess-Zumino action [l|. To avoid confusion we note the in this context we arc not talking about a 
Wess-Zumino term in the action of our (1+1) field theory that would have originated from integration of a closed 
3-form over a relevant 3-cycle. We rather interpret the potential in Eq. (f5"B"|) as the action of a free particle, represented 
by V e (with y playing the role of time), with an additional multi- valued term (ps, obtained via integration of a closed 
2-form over a relevant 2-cycle, as described by Eq. (|56|) . The quantity <£>b(ti) is the area enclosed by the loop n(y). 
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More precisely, values v?s(n) calculated with help of two paths x{Vi s ) an d ~x!{Vi s ) can differ by Anm with integer m 
if the paths belong to different equivalence classes. The notation <ps indicates that it is the Berry phase [j| associated 
with the loop n(y) on a sphere. Indeed, the parallel transport of a tangent vector on the unit sphere along the loop 
n(y) rotates the vector by the angle equal to the area enclosed by n(y). The Berry phase appears in the related 
quantum-mechanical phenomenon: the state of the spin with the projection s along the magnetic field acquires a 
phase factor e lS¥3B when the direction of the slowly changing magnetic field makes one turn of the loop n(y). 

Finally we note that the purely topological nature of driving, i.e., dF = 0, in the continuous field-theory limit of the 
considered model can be violated by a regularization, represented by a (finite) 2iV-dimcnsional system of N classical 
spins, arranged in a circular chain. Therefore, there is a question of how (and if at all possible) to regularize the second 
term in the expression for the driving force Eq. pip to preserve the topological nature of driving on the regularized 
level. This can be achieved by implementing an approach based on Eq. (|55[) . To define the topological term we note 
that in the thermodynamic limit N — > oo, the values of the neighboring spins are close, i.e., — rij\ <C 1 and, 

therefore, we can uniquely connect rij to rij + i with geodesic lines. This results in a piece-wise smooth loop n(y) with 
the desired topological potential uifB{n). We provide some mathematical details of the regularization in Appendix lAl 

B. Instantons and tubular neighborhoods 

This Subsection contains a preliminary discussion of our strategy in dealing with the problem of weak noise. The 
approach consists in reducing our original model to an effective stochastic model on a circle, where the new reduced 
variable is the re-parameterized Berry phase 9(x) £ S l : 

p B (n) = 27r(l-cos(0/2)). (57) 

The representation in Eq. (|57[) is possible and useful due to the instanton (optimal fluctuation) character of the 
current generation in this weak noise limit considered here. Then, a typical configuration is a closed loop on the 
sphere S 2 (also referred to, here and later on, as a "string") that is almost shrunk to a point performing a diffusive 
random walk over S 2 However, this typical diffusive meandering does not generate a current. Instead, the current 
is generated by rare events naturally occurring along the instanton trajectories illustrated in Fig. [1] The instanton 
trajectory is an optimal fluctuation as it has the highest probability (i.e., minimizes the Onsager-Machlup action) 
among the trajectories resulting in the transition. An instanton trajectory starts with a string, originally shrunk to 
a point, opening up into a "plane" circle configuration and then shrinking back to a point at the opposite end on the 
sphere. Such a fast and rare process generates a full cycle in the space of the re-parameterized Berry phase 9. 

The first step of our computational strategy, detailed in the next Subsection, consists in adopting the instanton 
approximation to evaluation of the current density J. The approximation means that the non-equilibrium current 
density J obtained from Eqs. (|41|) — (|43|) is concentrated in the narrow tubes near the configurations passed by the 
optimal fluctuation trajectories, i.e., the configurations x represented by "plane" circles embedded into S 2 . The 
concept of tubes where the current density is localized was used, for example, in Refs. [ill. [i~2l, [241. l25l, l39l. l43l. \5l\ . 

Let Mo C M be the subspace of these configurations. We refer to this as the space of instanton configurations. 
For example, all four configurations shown in Fig. [1] belong to the instanton space Mq. The instanton space Mq 
has a simple structure that allows parametrization for the instanton as well as for a small neighborhood C/ curr D Mo 
where the current is generated. First of all with a minimal abuse of notation we denote by dMo C Mq the set of 
configurations x G Mo with 9(x) = 0, i.e., constant loops x(y) = n, parameterized by their positions n £ S 2 on 
the target sphere. Obviously, dMo = S 2 . Since dMo C M is a compact manifold embedded into the configuration 
space, it has a standard small tubular neighborhood Uo D dMo, whose points (n, £) G Uo are parameterized by the 
position n on the sphere and the set £ of transverse variables. We further show that Mo \ dMo = SO(3) x (0, 2ir) is 
an embedded 4-dimcnsional non-compact manifold, thus having a small tubular neighborhood U D Mo \ dMo- This 
view suggests the following reparameterization of a "plane" circle x : S 1 — > S 2 . Let e 3 be a unit vector, orthogonal 
to the circle plane, with the direction determined by the loop orientation, e.g., by d y x(0). We denote by e\ a unit 
vector in the direction from the circle center to its origin a;(0), and set e2 = [e3, ei]. Then the oriented orthonormal 
basis set (ei, e 3 ) represents an element g € SO(3) of the orthogonal group. The fourth coordinate of x is its Berry 
phase 9(x) <G (0, 2tt). The points of the corresponding tubular neighborhood are parameterized by (g, 9, £) with £ 
representing the set of transverse variables. Obviously, f/ curr — Uq[JU covers the instanton space Mq. 

The aforementioned coordinate representation in the tubular neighborhood U has the following explicit represen- 
tation 



n(y) = sin(0/2) cos yd - sin(0/2) smye 2 + cos(6»/2)e 3 + °) e ^ 

i a 



(58) 
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where £ = = 1,2,...) is the transverse deviation expanded in the transverse modes with the components ipj(y) 
that generally is parametrically dependent on the Berry phase 8. For a regularized version j = 1, . . . , 2N — 4. We 
further note that since < 8 < 2ir in the region U, this region does not contain nontrivial 1-cycles. Therefore, the 
multi-valued potential Vwz restricted to U can be represented by a single-valued function, which can be expanded in 
the transverse variables as 

Vwz(x) = V (6)+W(6)(t<g)t)/2 Vo(0) = (tt/2)(-vcos<9 + 4ucos((9/2)), (59) 

where W(8) is represented by a symmetric matrix with the matrix elements Wij(0). 




FIG. 3: Illustration of the Vo(8) profile. Defined on the [0, 2ir] span the potential is multi- valued. The bars show the overlapping 
harmonic (blue) and WKB (green) domains. 

The most important for us is the ^-component of the force directed along the instanton trajectory 

F o (0) = -deVo(e) = (tt/2)(-« sin 8 + u ^2(1 - cos ffj) . (60) 

One observes that the effective potential Vq(6) has two stationary points, = and 8 = 0q where Fq(0) = Fq(6o) = 0. 
The unstable point 

= 2 arccos(u/w) (61) 

is characterized by 

fco = d e F {9 ) = {tt/2)v (1 - u 2 /v 2 ) > . (62) 

The stable point 8 = can be approached along two directions corresponding to different loop orientations. In terms of 
6 G (0, 27r) they correspond to 6 — > +0 and 6 — > 2tt — 0; in both cases dgF < 0. These two stationary points correspond 
to saddle-point {8 = 8q) and equilibrium (8 = 0) string configurations. The main instanton approximation, valid in 
the low-noise limit, means that only special instanton configurations, parameterized by 8 and possibly accounting for 
some small fluctuations (these from a small neighborhood U D Mq \ 8Mq), will be relevant to the discussion below. 

We will make some additional (to the basic instanton approximation) assumptions limiting the domain of validity 
but also adding a required extra tractability in exchange. In the following we will discuss separately: (a) The so- 
called WKB approximation that ignores the terms associated with second-order derivatives over 8 in comparison with 
the corresponding first- and zero-order terms. The WKB approximation is valid in the WKB domain where 8 is 
sufficiently far away from the equilibrium and saddle points: 8, 2tt — 6 \J n/k and \8 — 8q\ y/n/ko, respectively. 
The WKB region naturally splits into two sub-domains, < 8 < 8q and 8q < 8 < 27r. (b) The harmonic/ linear-force 
approximation is valid in a relatively small vicinity of the saddle-point, 1 6* — | 1, referred to as the harmonic 
region, where one can use the linear approximation for the force, Fq(8) « k${8 — do). Similar to the WKB treatment in 
quantum mechanics, the WKB and harmonic regions do overlap. See Fig. [3] for illustration. The concrete form of the 
instaton for the spin-chain model in the harmonic and WKB domains will be derived (and matched) in Section IV CI 
and Appendix [Bl Finally we note that, as will be demonstrated in Section IV C[ detailed analysis in the harmonic 
region of the stable point 6 = can be avoided. 
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Derivation of the Cramer function 



In this Subsection we apply the strategy, outlined in Subsection IIV C[ to the spin-chain model considered in the 
weak-noise limit also requiring that the topological currents are not too strong. As outlined in Subsection IIV C[ the 
vector potential A is characterized by the topological parameters given by Eq. (|4"T)l . As demonstrated in Appendix \X[ 
for the spin-chain model H±(M) = Z, i.e., there is only one topologically independent cycle s, which can be chosen to 
be restricted to the instanton space M . The cycle s for such a choice is the one illustrated in Fig. ((T|) and corresponds 
to just a change of 9 from to 2tt. The dual cycle a of codimension 1 can be chosen to be determined by the condition 
9(x) = 9o, so that s and a intersect at the saddle point. The topological current for our model is single-component 
and will be denoted by lu. The topological parameter is also single-component, 

Z = cxp ffiT 1 J^Ai{x)dx^j . (63) 

The four-step procedure for calculating S(u>) in the relevant for our application case dF = of the topological 
driving, was outlined in Subsection IIV CI Implementation of this procedure to the spin-chain model in the low- noise 
limit is facilitated by the following model-specific assumptions for the A, J, and p functions which can be verified 
directly once the solution, based on these assumptions, is found: 

(a) The vector potential A is essentially nonzero only within the sub-domain \9 — 9$\ < ^n/h®, which is contained 
in the harmonic region of the saddle point. In this domain A does not depend on the transverse variable C,, and 
its only substantially non-zero component, Ag, is along the 9 variable. 

(b) The eigenvalue A is exponentially small and can be neglected in Eq. (|49p everywhere except for a small vicinity 
of the equilibrium point 9 = 0. In this sub-domain A is so small that it can be totally neglected. Moreover in 
the entire WKB region A is still small enough, so that the nonlinear term A 2 in Eq. (|4"5|) can also be neglected. 

(c) The main contribution to the integral for S(u>) in Eq. (|52j) comes from the | ^ — 6*o | ^ y/t/fco domain. 

(d) The current density J is concentrated in a small tubular neighborhood of the instanton space Mo, whose 
transverse size scales ~ ^fn with the noise value. In the WKB and saddle-point harmonic regions the longitudinal 
component of J has a Gaussian dependence on £ 

Jg = J Q (9)e-^ laW ^\ (64) 
In the harmonic region of the saddle point this is the only non-zero component of J, and a(9) ~ W(9q). 

(e) The density distribution p is concentrated near the equilibrium 9 = 0, where rare events generating the current 
can be neglected: 

p(x) ^ p e- 2 ^ 1Vwz(x} . (65) 

The features of A, J, and p, listed above, are generic for a topologically driven system in the low-noise limit. 5*0(3) 
symmetry is an important special feature of our model. Therefore, the critical points are represented by isolated orbits 
of 50(3) rather than isolated points. Discussing "the" equilibrium 9 = and "the" saddle 9 = 9q points in our model, 
what we actually mean is that we have two isolated orbits of 5*0(3) which represent the equilibrium and saddle points, 
respectively. These orbits are given by 50(3)/50(2) = 5 2 and 50(3), respectively. The symmetry gives rise to zero 
modes, which substantially complicates a straightforward path-integral calculation in the instanton approximation, 
especially due to the different numbers, 2 and 3, of zero modes at the equilibrium 9 = an transition 9 = 9q 
configurations, respectively. 

In our approach, the zero modes are allowed for automatically in a very simple way in the form of the volumes of 
the relevant orbits 5 2 and 50(3) for the equilibrium and saddle configurations, respectively. 6 



Our spin-chain model is degenerate and thus special, in what concerns the number of stable points (only one) and number of zero modes 
at the stable and unstable fixed points. In a generic system the number of the equilibrium and saddle points can be arbitrary. Moreover, 
the instanton manifold is 1-dimensional. An instanton trajectory from the manifold starts at a stable point ascends to a saddle point 
and consequently descends following an unstable direction to another stable point. Therefore, Mq in general can be viewed as a graph 
of instanton transitions. Using the four-step strategy, outlined in Subsection IIV Cl the problem of finding the Cramer function of u> can 
be reduced to a Markov chain model on the graph that represents the instanton space Mq. This general approach, allowing to reduce a 
non-equilibrium field-theory problem in the weak noise limit to a Markov chain model with states associated with stable fixed points of 
the classical dynamics, will be addressed in a separate publication. 
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Properties (a)-(e) listed above allow us to ease computation of the Cramer function S(u>) essentially and in particular 
bypass complications associated with detailed resolution of the WKB domain. Therefore, in this Subsection we will 
rely on these properties, presenting justification details related to WKB calculations, as well as some technical details 
on the determinant calculations, in Appendix IBl 7 

The material in the remaining part of the Subsection is split into paragraphs according to the four-step strategy 
described in Subsection IIV CI 



1. Step (i): Identifying the vector potential 

In this Subsection we implement step (i). In the harmonic region \9 — 6q\ <C 1, according to property (a), the 
vector potential A is described by the only non-zero component Ag that depends on 9 only. We seek a solution in 
the representation of Eq. flU)) with p-{9) depending on 9 only. According to property (b) we set A = in Eq. P5j) . 
Making use of Eq. (|46[) we arrive at the following homogeneous linear equation: 



(«/2)$p_ + F (9)d e p- = , (66) 
whose general solution can be expressed in terms of the error function erf(z) = 2tt~ 1 / 2 Jq dsexp(— s 2 ): 

P-(9)/p-(9 ) = 1 + cerf (y/hfa(6 - o j) ■ (67) 
This results in the following expression for the ^-component of the vector potential: 

a u I i cexp(-fc (6>-6>o)7 K ) 

A g = -^Anko/it (68) 

1 + cerf (yV^(0-0 o )J 

To determine the constant c we recast Eq. (|63j) as Z = exp^J^ d9Agj re), which results in 

c=(l-Z)/(l + Z). (69) 

Eq. (|68p shows that A peaks near 9q and decays in a Gaussian fashion as 9 moves away from 9$. The solution 
presented in Appendix IB 11 shows that the rapid decay continues with 9 moving even further away from 9q inside the 
WKB domain. 



2. Steps (ii) and (Hi): Identifying the current- density and density distributions 

The current density J does not have any special structure in the harmonic region in contrast to the distributions 
A and p. Therefore, on the step (ii), the solution of Eqs. ([50)1 for J near the saddle point can be extrapolated from 
the WKB region (see Appendix IB 2ft . 

Step (iii) starts with solving the first equation in Eqs. (j5Tj) . This is an easy task, since Eq. (|5T|) . being restricted to 
any 1-dimensional subspace, becomes a linear first-order ordinary differential equation with a right hand side. On the 
cycle s that belongs to the instanton space, the density p is represented by a function g(9), and the linear equation 
adopts a form 

( K /2)dgg(9)-(F (9)-Ag(9)) g (9) = J e (0), g(0) = g(2n) = Po (70) 

The 1-dimensional linear boundary problem, defined for < 9 < 2tt and represented by Eq. (|70p . can be solved in a 
standard way [39[ by using a representation 

g g 

g{9) = q{9) cxp (-2 K ~ l V{9)) , V{9) = V {9) + J Ag{9')d9' = J (Ag(9 r ) - F (9'))d9' (71) 



7 Note also that the shape of the current density distribution J derived in Appendix [B] provides with some useful information on the 
processes that generate the current. 
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Substituting Eq. (|7T|) into Eq. ([70)) one arrives at 

d e q = -(2/ K )J e 2K ~ 1 v( e \ q(0) = p O) g(2^)=p expf2«- 1 / d0(4 e (i 
where we have introduced the barriers 



(i 



S+ = -2k- 1 J d6 F {6) , S-=2k~ 1 J d6F (6). (73) 

O 

Close to the saddle-point configuration, q depends only on and the only substantially nonzero component of the 
current density J is 

J e (x) = Joe-"" 1 * 8 * , (74) 



where a = W(9o) and Jo is constant on the length scale y 7 n/ko, according to the property (d). 
Integration of Eq. (|72|) results in 



q{9) =p - -J e s + / drexp I -k (r - 9 ) 2 /k+(2/k) I dr 1 A 6 {t') I , (75) 



H 



o \ o 

which translates with the help of Eq. (|68f into 



q{9) = PQ - ^/{^Joe^ — + ^/{^ko)J e s + — (1 — , (76) 



where c = (1 — Z)/(l + Z). One observes from Eq. ([76]) that g(6*) is localized at \9 — 9q\ < -Jn/ko, thus justifying the 
approximations used so far to evaluate Ag and Jg. Setting 9 = 2-k in Eq. f|T6[) and applying the boundary condition 
from Eq. (|72|). one expresses Jq in terms of Z and pq: 



Jo = poV^o/(47r)(Z- 1 e- s + - Ze- S -) . (77) 

To complete the major step (iii) we apply the normalization condition, given by the second relation in Eq. (j5ip , to 
find the normalization constant po. We consider the neighborhood Uq represented by the most probable configurations 
that are close to constant loops. Stated differently, the density given by Eq. (|65[) is concentrated in the harmonic 
region of the stable point and hence adopts a form: 

p{x) = p{n,i) = poe-^ 1 ^ 3 , (78) 

where n G S 2 is the center of a small loop, and the matrix 7 determines a harmonic expansion of Vwz(x) hi terms 
of nonzero modes Then the normalization condition becomes 

f dtt(n)Jl>Sp(n,t) = l, (79) 

thus leading to the following value of po (see Appendix IB 31 for more details): 



po = (l/4) v / d^(7r^)- JV K 1 - Ar . (80) 

Note that due to the symmetry of the model dp,(n) is an 50(3)-invariant measure on the sphere, defined up to a 
multiplicative factor. This factor has been identified in a standard way by considering the two zero modes on the 
sphere (see Appendix IB 31 for some detail). The determinant det7 includes (2JV — 2) positive eigenvalues of the 
discretized operator 7. 
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3. Step (iv): Finding the Cramer function for the topological current 



Step (iv) starts with evaluating the integral (|52[) for the Cramer function iS in terms of Z. According to the 
property (d) the integral acquires its major contribution from the vicinity of the saddle-point configuration. Indeed, 
the integrand decays as oc exp((2/re)Vo(0)) with 9 deviating from 9 inside the WKB domains, where Eqs. (168p . 
(|69|) . (|7Tj) and (|76|) can be used. The dependence of p on the transverse variables in the relevant region is given by 
p(x) = g(9) exp(— K _1 cr^ ® £) with a = W(9 ), which follows from Eq. (fSTj) applied in the transverse direction as well 
as from the asymptotic absence of the transverse components of A and J in the relevant region. 

Thus, we obtain the following expression for the Cramer function: 



2jv 



S = (2k)- 1 f dfig (g) fvCe-^ 1 ^ f d9A 2 e (9)g(9). (81) 

JSO(3) J J 



The integral over 9 can be calculated similarly to that in Eq. ([75]) whereas the calculation of the other integrals is 
presented in Appendix IB 31 



S = d Q N N {irK) N - 2 /Vdct<7 (-J Q \nZ + ^Kko/irpoil - Z)(Zer s - - e" s+ )/(2Z)) . (82) 
The integral over the space 50(3) of the zero modes is computed explicitly: 

O = / <Wff) = 2tt 2 (1 + u 2 /v 2 )^l - u 2 /v 2 , (83) 

JSO(3) 

where the invariant measure dpg a (g) in 50(3), taking into account three zero modes produced by infinitesimal changes 
of g, requires careful evaluation of the multiplicative factor. This is performed in a standard way by considering the 
scalar products of these modes (see Appendix IB 31 for some detail). 

To conclude the step (iv) we need to relate the topological current lo to the topological parameter Z by applying 
the second relation in Eq. ((52)) . Let us reiterate that we have chosen the cycle a by the condition 9(x) — 9q, so that 
we have for the current 



which results in 



cj = u>*[a]= / ,7 = / dxJ e (x) , (84) 

Ja Je(x)=e a 

w = Jotfo J VCe-^ 1 *^. (85) 

After calculating the Gaussian integral in Eq. (|85p over (2N — 4) modes £ with positive eigenvalues at the saddle 
point (cf. Eq. (|8"Tj) ). one derives 

u = Jotf Q N N (irK) N - 2 /VdctcT. (86) 
We combine Eqs. (|80|) . (|77|) and (|86|) to get the topological current in terms of Z: 



lo = iV(4tt 2 k) v/det 7/ det ct V 'nk / '(4 7 r)(Z- 1 e" s + - Z e - S -) . (87) 

This equation has a form 

w = - (88) 

where the newly introduced quantities 

k± = tV(4tt 2 k) y/det 7/ det crv/Kfco/(47r)e" s± (89) 

can be interpreted as the transition rates in the auxiliary Markov chain model discussed in detail in the next Subsection 
IV C 41 The stationary current can be obtained by setting Z = 1. The final expression of Z via lo follows from Eq. 
(the choice of the root is related to the model equivalence discussed in Appendix ICj) : 
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4- Final expressions and two-channel single-state Markov chain 

We complete the derivation of the Cramer function by providing the relation to the parameters of the spin-chain 
model v and u. The ratio of the determinants can be calculated in the limit N ^> 1 (sec Appendix IB 31 for details): 

v/det7/dctcr = ■k~ 2 v 4 u' 1 (v 4 - it 4 )" 1 sm(iru/v) . (91) 
The final expression for the Cramer function is 

S = -wln-¥ — -± y/u 2 +4k + k_ + n+ + k_ , (92) 

ZK- 

where k± are expressed through v and u with the help of Eqs. (|62p. ([55)1 . (|89p. and (|9ip. whereas the barriers are 

S± =2tt(vt<)7M- (93) 




FIG. 4: Two-channel single- state Markov chain. 

Now we can demonstrate the equivalence of the spin-chain model to the simple Markov chain model, shown in Fig. 
[Hand described in Appendix [Cl and justify that A can be neglected in the above solution. Using the properties of the 
distributions J(x) and A(x), we can rewrite the first representation of the Cramer function in Eq. (|53[) as 

S = -uj\nZ - A. (94) 

The eigenvalue A can be calculated by comparing the two representations of the Cramer function: 

A = Z~ 1 n + + Zk_ — (k + + K-) , (95) 

or in terms of uj, 

A = \J <J 2 + 4k + k_ — k + — k_ . (96) 

The equivalence of the spin-chain model at moderate topological currents to the two-channel single-state Markov 
chain model is, in particular, obvious from the forms of Eqs. ([9"0")l . and ([9"5)l . 

We see that A can indeed be neglected if the total current is not exponentially larger than the equilibrium current. 

VI. SUMMARY AND CONCLUSIONS 

Let us briefly summarize main steps, results of the general approach and the model discussed in the manuscript. 

We have developed here a topological picture of generated stochastic currents, where each current component is 
associated with a topologically nontrivial 1-cycle in the system configuration space, so that the current u> £ Hi(M; R) 
resides in the first homology of the configuration space with real coefficients. The current, defined in a topological 
way as the set of number of rotations over the independent cycles per unit time, is related to the current density J via 
the Poincare duality H±(M;M.) = H m ~ 1 (M;M), where the cohomology is considered in the dc Rahm representation. 
( See. e.g., [32j for review of the de Rahm cohomology.) By considering the current density J as a closed form of rank 
(m— 1), the current is viewed u> = [J] as its cohomology class. This allows the Cramer function <S(<*>) of the long-time 
current distribution to be calculated by applying the variational principle to the current-density functional S(J,p). 
The explicit form of S(J, p) together with the variational principle and an intuitive path-integral base derivation is 
also presented in the manuscript. 

We further focused on the case of topological driving, when the driving force is locally potential, dF = 0. Even 
though the system may not be represented by a potential function due to presence of topologically nontrivial cycles, 
it still allows description in terms of a multi-valued potential F = —dVwz, referred to as a Wess-Zumino potential, 
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due to it close resemblance with multi- valued Wess-Zumino actions well-known in quantum field theory. By applying 
the variational principle to the calculation of the Cramer function S(u)) we derived a system of equations for p, «7, 
and A, the latter being an auxiliary curvature-free (dA = 0) gauge field. We also developed a procedure of solving 
the aforementioned equations step-by-step, resulting in the well defined algorithm for calculating the Cramer function 
iS(u>) of the topological current. 

To illustrate the general approach we considered a circular spin-chain stochastic model with topological driving 
that constitutes a regularized version of a (1 + 1) stochastic field theory. The model represents Langevin dynamics of 
an elastic string evolving over a two-dimensional sphere S 2 . Despite of its high intrinsic dimensionality the problem 
appears tractable as it has only one independent topologically nontrivial cycle, H\(M;M.) = R, which reflects the 
well-known result on the homology of the loop spaces of spheres. In the low-noise limit we have solved aforementioned 
variational equations for p, J, and A explicitly. The solvability became possible via proper use of the rare-event, 
instanton analysis of the current generation. Moreover, the actual calculations are streamlined even further, as in 
fact we do not perform the path-integral calculations around the instanton solutions explicitly, but rather use the 
instanton scenario to build an asymptotically exact ansatz for (p, J, A) . This useful technical trick allows us to solve 
the equations (in the proper low noise and large deviation limit) analytically. In particular, we have found that 
the current density J, generating the topological current u>, is concentrated in a small tubular neighborhood of the 
finite-dimensional space Mo C M of instanton configuration, and it shows a Gaussian dependence on the set £ of the 
transverse variables. Technic ally , thi s appro ach can be viewed as an extension of the equilibrium techniques based 
on current tubes [Til . 

m, m, m, m, m m to the non-equilibrium problem of evaluating the Cramer functions of 

generated stochastic currents. It is instructive to note that, although the relevant distributions are concentrated in a 
small tubular neighborhood of the whole instanton space M , the calculation of S(<u) requires careful consideration 
of a small neighborhood of the saddle point (transition state) only. This reflects the fact that, during a rare event 
that contributes to the current generation, the system spends most of the time around the saddle point, being thrown 
there by a strong fluctuation of noise and waiting for a small extra kick that starts its unavoidable fall-down back to 
the steady state. Most of the time the spin-chain is involved in a "boring" , i.e. typical, diffusive meandering along 
the sphere, thus waiting (almost forever) for the next instanton jump/transition. This scenario allows us to reduce 
the complex stochastic field theory to a simple Markov chain model. 

This asymptotic (low-noise, instanton) reduction of a continuous stochastic system to a simple Markov chain model 
in the case of moderate non-equilibrium currents is quite a general result. In the case of a purely potential force 
F = —dV it is connected to the topological properties of the potential and closely related to the Morse theory [34[. A 
connection between the Morse theory and instantons was established in (ol| , where a super-symmetric imaginary- time 
quantum mechanics has been introduced with the effective potential V e s = (dV) 2 . In the semiclassical limit h — > 
the instanton approach has been implemented to evaluate the effects of tunneling between the metastable states, thus 
allowing us to distinguish between the true zero modes, which contribute to the de Rahm cohomology (32l . [44j and 
hence describe the topological invariants of the configuration space, and just soft modes with exponentially vanishing 
eigenvalues in the h — > limit. This approach is known as the Morse- Witten (MW) theory [61(. It is well known 
[37j that a simple gauge transformation turns the Fokkcr-Planck (FP) operator into the Schrodinger operator in the 
potential V e g, with h 2 representing the temperature in the FP theory. Within this equivalence the instantons have 
the same shape, although in the FP picture they play the role of optimal fluctuations that minimize the Onsagcr- 
Machlup action for rare stochastic transitions between the metastable states. The eigenvalues of the soft modes, 
describing tunneling in the MW theory, attain a true physical meaning in the FP picture, since they represent the 
slowest relaxation rates, which are due to rare over-the-barrier transitions. Such an approach has been utilized by 
Tanase-Nicola and Kurchan [57| who studied a supper-symmetric FP theory that extends the standard FP dynamics 
in the same way as the super-symmetric quantum mechanics extends standard quantum mechanics. 

At this point we note that the standard Morse theory describes the case of a purely potential force, whereas for our 
spin-chain model the force dF = is only locally potential. The topological counterpart of stochastic dynamics in 
this situation is known as the Morse-Novikov (MN) theory that has demonstrated its capabilities to study topological 
properties of the underlying configuration space [47] . It is worth noting that in the continuous limit the multi- valued 
potential Vwz(n) of our spin-chain describes the multi- valued action of a free particle moving along the sphere S 2 
in a magnetic field of a constant curvature. If the particle motion is affected additionally to the magnetic field by 
a non-zero pot ential, the problem, which can be mapped onto the Kirchhoff equations, can be handled using the 
MN theory [471 ] . In a generic situation the instanton space Mq is one-dimensional and is represented by the unstable 
spaces of the isolated critical points of Morse index 1 (one unstable mode). The spin-chain model, considered in this 
manuscript, whose continuous limit corresponds to the Kirchhoff problem with the zero potential, is degenerate due 
to its intrinsic 5*0(3) symmetry, restored when the potential vanishes. In terms of the Morse theory this degeneracy 
means that the critical points, which arc isolated in the case of standard Morse (including MN) theory are replaced 
by isolated critical manifolds associated with the orbits of the underlying symmetry group, 50(3) in our case. Such a 
situation, in the potential force case is described by the so-called equivariant Morse-Bott (MB) theory Therefore, 
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the topological counterpart of the low-noise stochastic dynamics for our spin-chain model can be referred to as a case 
of the equivariant Morse-Bott-Novikov (MBN) theory. Stated in the stochastic dynamics terms, the instanton space 
is represented by a 1-dimcnsional family of the SO(3) orbits, actually 50(3) itself, almost everywhere. Then only 
the stable points (constant loops) are represented by an orbit from S 2 = SO(3)/SO(2). However, and in spite of this 
degeneracy, we have just showed that this collapse of the orbits at the stable critical points is not an obstruction for the 
general Cramer function calculations. Let us also emphasize that even though traditionally the standard (equivariant) 
Morse theory was developed primarily for the case of a potential force, we found out that in fact it can be efficiently 
used in the stochastic dynamics aspect in the case of moderate topological driving, when the Morse function is multi- 
valued. As we just demonstrated on the example of the non-equilibrium spin-chain, moderate character of the driving 
does not change the topological structure of the critical and instanton manifolds in comparison with the pure potential 
case. 

In a generic non-equivariant case with moderate topological driving one could expect a full reduction of the low- 
noise stochastic dynamics to a Markov chain process on a graph whose nodes and links represent the critical points 
of Morse index and the unstable manifolds of the critical points of Morse index one, respectively. This is really the 
case when the potential function satisfies the Morse-Smale (MS) condition 34[, which would be a generic situation. If 
the MS condition is not satisfied (which is some kind of degeneracy) the stochastic dynamics of the current generation 
can show some additional interesting features, which are yet to be analyzed. In the case dF ^ of intrinsically non- 
potential force, whose topological counterpart is represented by Morse decompositions in the Conley index theory [13], 
stochastic dynamics can be extremely complicated, with the critical spaces, represented by neither isolated points nor 
even isolated manifolds but rather some closed sets, possibly of fractal nature. Even in the simplest non-equivariant 
case of isolated critical points, the situation is much more complicated (than in the equivariant case), yet apparently 
treatable. The difficulty is that transitions between different isolated states of the effective Markov chain become 
direction-dependent, since as opposed to the potential and locally-potential cases (with isolated critical points) the 
instanton trajectories that correspond to climbing the barrier and falling down the barrier are different. The problem 
of computing the pre-exponential factors in this case also becomes much more involved and less universal. 

Let us also stress that in the v = case the Onsager-Machlup action of the spin-chain model in the thermodynamic 
field theoretic limit, N — > oo, reproduces the action of a (1 + 1) sigma-model with the topological term. However, 
in the regularized model (finite N) the relative fluctuations in the v = case are strong, so that rij+i is typically 
not close to rij. The corresponding string on the sphere is not continuous, which can be viewed as the main reason 
for complexity of the sigma-model, considered as a field theory. The clastic term proportional to v suppresses the 
relative fluctuations in such a dramatic way that statistically rij+i and rij are always close to each other. Therefore, 
the configurations can be viewed as continuous loops in S 2 even for a finite large N, which creates a topologically 
nontrivial cycle. This interesting peculiarity of the problem gives rise to generation of the single component stochastic 
current. We also note that since in the v ^ case and TV — > oo the result is finite, the elastic term in the case of small 
v can be interpreted as a regularization of the sigma-model based stochastic field theory suppressing the short-range 
fluctuations and eliminating divergences. 

The methodology for the LDP of empirical currents, developed in this manuscript, is closely related to the FP 
approach to MW theory developed by Tanase- Nicola and Kurchan [53] • In a way, this methodology can be interpreted 
as an extension of the approach of (57| to the locally potential driving case dF = 0, implemented in low-dimensional 
0- and 1-fermion sectors that correspond to the density and current density distributions, respectively, and applied to 
large-time statistics of empirical currents. Thus, we believe that our weak noise analysis of currents in the problems 
with topological driving can be restated in super-symmetric terms of [13] • Testing this conjecture remains a future 
challenge. It would also be interesting to extend the above relation to higher dimensions, e.g., by treating diffusion 
in n-fcrmion sector as noisy dynamics of n-dimcnsional chains that produces higher dimensional currents, residing in 
iJ„ + i(M; R). Based on the presented in the manuscript low-dimensional results, it appears that the low-noise limit 
of the Cramer functions of the empirical currents (including the higher-dimensional counterparts) contains detailed 
information on the Morse decomposition of the underlying space M. 

Finally, we have considered the simplest stochastic (1 + 1) sigma model with the target space represented by S 2 = 
CP 1 . A simple topological computation demonstrates that a generalization based on the target spaces CP n with n > 1 
still provide a single component topological current. The simplest generalization resulting in the multi-component 
topological currents u) requires the target space to be a complex flag space [H, H3]. The Cramer function S(uj) can 
be computed in such cases straightforwardly and without any significant complications, by using the methodology 
explained above. 
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APPENDIX A: TOPOLOGICAL STRUCTURE OF THE CONFIGURATION SPACE 

In this Appendix we discuss mathematical formulations for the spin chain regularization introduced in Section IV Al 
in a somewhat lighter, more physical terms. 

Let us first introduce a simple and useful representation of the first homotopy and homology groups of the relevant 
configuration space Map(5 1 ,S' 2 ) of maps from the base space 5* 1 to S 2 . To regularize our statistical field theory we 
should consider the spaces of maps S 1 — > S 2 , with various degrees of smoothness. Ultimately, we are interested in 
random maps (non-smooth) whose topological properties can be approximated using well-behaved (continuous) maps. 

The spaces of smooth and piece- wise smooth maps from S 1 to S 2 are denoted by Map°°(S' 1 , S 2 ) and Map^(5 1 , S 2 ), 
respectively. The topology of these map spaces is defined in a standard way. To analyze proper discretizations of 
the spaces of smooth maps, we represent the circle S 1 by a cyclic lattice j = 0, . . . , N — 1, approximate a map n(y) 
by a set of points rij, and introduce e = 2ir/N. Finally, we define a set of approximations Ln i£o S 2 for the space 



close to each other. For not too large e , e.g., for e < 1/3, we can define a continuous map L N _ eo S 2 — ► Map(S' 1 , S 2 ) 
by connecting the neighboring points rij, rij + i with geodesic lines. 

This map generates homomorphisms between the homotopy, homology, and cohomology groups, respectively, for all 
the three spaces of maps. It is possible to show that if an approximation is accurate enough, the relevant topological 
properties of the approximations Lpf t£o S 2 are identical to these of the original space Map(5' 1 , S 2 ). 

To identify a very simple and intuitive picture of the topological current generation in our field theory, we define 
continuous maps 9 : M -> S 1 for M = Map°°(S 1 , S 2 ), M = Map^°(5 1 , S 2 ), and M = L N ^ Eo S 2 by associating with 
any loop from M its Berry phase defined as the holonomy along x £ M, also understood as x : S 1 — > S 2 . 

Introducing a map from e tVB/2 £ U(l) to 9 £ [0, 2tt] S 1 and combining it with the maps defined above, we arrive 
at the continuous maps 6 : M -► S 1 for M = Map°°(S\ S" 2 ), M = Map~(S rl , S* 2 ), and M = L N ^ Eo S 2 . In all three 
cases the morphisms ni(M) tti(S 1 ) = Z and H\{M) — » Hi(S 1 ) = Z generated in the homotopy tti and homology 
Hi groups, respectively, are isomorphisms. All this implies that the current generation can be observed by monitoring 
the reduced variable 9(x) £ S 1 and counting the windings around the circle. 

APPENDIX B: DERIVATION OF THE CRAMER FUNCTION FOR THE SPIN-CHAIN MODEL: 

DETAILS 

In this Appendix we present some details of the Cramer function <S(u>) derivation for the spin-chain model. 



Since the WKB region does not contain topologically nontrivial 1-cycles, Eq. (JBTJ) can be recast in the following form: 




1. Vector potential in the WKB region 



In the WKB region we set A = and neglect the nonlinear term in Eq. (|49|) , which turns it into 

(n/2)S A + F ■ A = 0, dA = 0. 



(Bl) 



=0, A = dip 



(B2) 



We seek the solution of this equation in the following form 



ijj-(x) = ip(x)e 2K ly »-^ x \ ifj( x ) = V>o(#)e 



— 1 



"W(C®0 c^(x) = 0, 



(B3) 
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or, equivalcntly, 

4>-(x) = ^ (6))e 2K_ll/o(e) - K_1<T - (e)(c ® c) , a-(0) = a(6) - W{9). (B4) 

Wc further substitute tJj(x) given by Eq. (|B3[) in the equation Cijj(x) = 0, retain only the leading terms in k, and 
keep in mind that typically |C| ~ \f^- Overall this results in the following system of equations: 

- d 8 (F ip ) - V'oTr(cr - W) = 0, -F a Vga = 2cr 2 - Wa - aW. (B5) 

Combining the two Eqs. (|B5[) one derives 

9e(FoVo/Vdet a) = 0, -F V 9 tT = 2cr 2 - Wcr - cjW. (B6) 

The second equation in Eqs. (|B6j) is nonlinear, and thus intractable, if the matrix is large. However, one can still 
show that this equation does have a solution with the following important property, a(8) = W{9) for 9 = 0, 9 , 2n. 
To demonstrate this property one refers to a dynamical equations with respect to some proper time r 

6 = F Q (6), a = -2a 1 + aW{6) + W(9)a, 6 = d T 9, a = V T a. (B7) 

This dynamical system has three critical points (9,W(9)) with 6 = 0, 9q, 2tt, where 9 = 0, 2ir correspond to stable 
critical points, while 9 = 8q describes an unstable critical point (deFo(8o) > 0). Therefore, finding a solution with 
the desired properties is equivalent to identifying two solutions of Eq. (|B7[) which start at the unstable point and 
reach the two stable points in infinite time. Generally such trajectories do exist. For our spin-chain model it can be 
demonstrated in a straightforward manner, e.g., by using the angular-momentum representation for the transverse 
modes tpj(y;9) entering Eq. ([58]) . 

The first equation in Eqs. (|B6[) can be easily solved: 

M0) = Cy/dct<r(0)/F o (O), (B8) 

thus immediately providing expressions for Ag, via Eqs. (|B2j) and (|B4j) . Therefore, the WKB solution for Ag that 
matches with the harmonic solution (|68[) within the domains where both apply, has the following form in the two 
distinct sub-domains: 

0<#<#o: Ag = -^ K k deta/(irdetW(9 ))(l/Z - i) e 2 «" 1 ( v oW-v (eo)) e -«- 1 <r-(e)(C®C) ) (B9 ) 
0o<0<2n: Ag = -v / Kfcodeta/(7rdetH/(6»o))(l-^)e 2 ^ 1(yo(e) - Vb(eo)) e- K_1 ' T - (e)(c ® c) . (B10) 
The expressions for the transverse components Aj of the vector potential can also be derived straightforwardly. 



2. Current tubes in the WKB region 

In this Subsection we determine the current density distribution J. To achieve this goal we solve Eq. (|50p assuming 
that J is concentrated in a small neighborhood of the instanton space Mq. 

In the WKB region we can set A = in the first equation in Eqs. ((50|) . and since the region does not contain 
topologically nontrivial 1-cycles, the first equation is equivalent to J = {n/2)dtp — Ftp for some function <p(x). 
Substituting this representation into the second equation of Eqs. ([50)1 one arrives at the following system of equations 
for J 

&p(x)=0, J = (n/2)dip — Fip. (Bll) 

In the low-noise limit the current is generated by rare events, and during the transition the system experiences small 
Gaussian fluctuations around the instanton trajectories. Therefore, the current density J is concentrated in a small 
tubular neighborhood U cum - D Mq of the instanton manifold Mq where its longitudinal current component should have 
a Gaussian dependence on the transverse variables £• To verify this property we seek the solution of Eq. (|B11[) in a 
form 

cp(x) = <p(0, C) - M0)e- K ~ la{9){Cm) - (B12) 

We substitute the ansatz of Eq. (|B12|) into Eq. (|B11[) and apply the same strategy as used earlier in Section IB II in 
the context of deriving Eqs. (|B5|) and (|B6[) . This results in the following system of equations 



d e (F <p /y/deta) = 0, -F V e cr = 2cr 2 - Wcr - aW. 



(B13) 
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Since the second relation in Eqs. (|B6[) and (|B13[) are identical, using the same notation a for the variances in Eqs. 
(|B3|) and (|B12[) is perfectly legitimate. 

Solving the first equation in Eq. (|B13j) . further substituting the obvious solution into Eq. (|B11|) . and keeping the 
leading terms in ^[k in the way detailed in Section IB 11 we arrive at the following expression for the current density 
J in the WKB region: 

Je(0,C) = Jo Vdet a{9)/ det W(9 ) e -^ 1 ^^\ 

Ji(e,C) = JoV^^)/^w(e )(F (e)y\a lJ (e)^w lJ (e))Qe~^ 1 ^ 9 ^<l (bu) 

Careful examination of the second equation in Eqs. (|B13|) in the harmonic region of the saddle point and allowing 
for the properties of the solution described at the end of Section |B II show that <J-{0) = a(Q) — W(9) tends to zero, 
when 9 approaches 9q, as ~ (0 — Oq) 2 . Therefore, the transverse components Jj vanish in the harmonic region. This 
confirms the assumption we have made in Section (|V C[) . referred to there as the property (d). 

At this point we note that the current distribution J, given by Eq. (|B14[) derived for the WKB region, also extends 
nicely into the harmonic region of the saddle point 6 = 9o, thus suggesting that J in this region is also described by 
Eq. (|B14|) . This can be verified directly. The reason why the WKB solution easily extrapolates into the harmonic 
region is related to the asymptotically longitudinal nature of the vector potential (A, = 0) in the harmonic region of 
the saddle point. 



3. Computation of the relevant determinants 

In this Section we present some details of the functional integral calculation in Eqs. ([75]) . (fSTj) . and In Eq. (17^)) 
the integration is performed over the deviations from the stable constant loop, whereas both in Eq. (|8 1 [) and Eq. ([84]) 
the integral runs over the transverse deviations from the saddle-point loop. In the following calculation, as in the main 
text, both the density p and the current density J arc assumed regularized on a lattice of N spins, in accordance with 
discussion of Appendix [SJ Therefore, the functional integrals should be represented by finite-dimensional integrals 
over N positions rij on the unit sphere. Since we will see that only long-wavelength deviations from the constant and 
saddle-point loops contribute to the ratio of the two integrals of interest, we will simply execute the limit N — ► oo in 
all the intermediate expressions where it exists. 

The integral in Eq. (|79p is evaluated over the two-parametric deviations 6n a (y) (with a = 1,2) from the constant 
loop. The potential accounting for configurations around the constant loop is harmonic: 



V WZ « (l/2)(*n,7<*n) , 7 = 2tt ( ) , (B15) 



where the scalar product is conventionally defined as 

r-2n 



«,»?)= r ^££foto«(y). (Bi6) 

Jo 2tt^ 

The operator 7 can be diagonalized in the space of Fourier harmonics e lqv with q = 0, ±1, ±2, . . . , ±(N — l)/2. For 

the sake of convenience, and since it does affect the N — > 00 limit, we consider an odd N. The eigenvalues and the 
corresponding normalized eigenvectors of 7 are 

j q = 2n(vq 2 + uq) and j q = 2n(vq 2 — uq) , (B17) 

E q = (1/V2)e^ ( \ ) and E q = (1/V2)e^ (]). (B18) 



There are two zero eigenmodes in this set: .Eo and Eq. All other eigenvalues are positive in the considered case v > u. 
One arrives at the following expansion in terms of the eigenvectors: 



fzN) J2( C * E * + ■ ( B19 ) 



Since 5n is real and E* = E_ q , the coefficients are related as c_„ = c*, and the transformation from the set of 



27V spin vector components {(8n±(y), Sri2(y))} with y = 0, £,..., 2ir — e (where e = 2it/N) to the set of N Fourier 
harmonic components {(Rec g , Imc g )} with q = 0, ±1, . . . has the Jacobian equal to unity. The zero mode coordinates 



2G 



{Reco,Imco} are related to a uniform shift 5n = const along the sphere as cq = VN(Sni + 18712). Thus, the integral 
over all deviations 8n from the constant loop consists of a zero- mode factor AirN and 2 (AT — 1) integrals over positive 
modes £, finally giving 

J dn{n) J D£exp(-KT 1 ($,'yt)) = 4n N N N k 1 *- 1 / ^dctj , (B20) 

where the determinant det7 includes 2 (TV — 1) positive eigenvalues of 7. 

The integral in Eqs. (|5Tj) and over transverse deviations C from the saddle-point loop no(y) and its SO (3) 
rotations can be calculated in a similar fashion. All small deviations, including the longitudinal ones and rotations, 
can be decomposed into meridional and zonal components: 

Sn(y) = S ni (y)§(y) + 6n 2 (y)j>(y) , (B21) 
0(y) = ei cos(6*o/2) cosy — e 2 cos(6*o/2) siny — sin(6*o/2) , cj)(y) = e\ siny + e 2 cosy, (B22) 

with the unit vectors introduced before Eq. (|58[) . The potential expansion of Eq. (|59[) can be represented in terms 
of the scalar product defined in Eq. (|B16[) as 

V wz (n + 5n) - V (8 ) = (l/2)(«Jn, <x6n) , a = 2tt ( " U " /v ^ ^ ) . (B23) 

Note that in this Appendix a acts in the space of all deviations including the negative mode 8ri\(y) = const = 56/2, 
whereas in the main text a is restricted to the subspace of transverse deviations corresponding to positive modes. 
The eigenvalues and the corresponding normalized eigenvectors of a can be found similarly to those of 7: 



tr,,± = 2tt [vq 2 - (v/2) (1 - u 2 /v 2 ) ± \] q 2 u 2 + (t> 2 /4) (1 - u 2 /v 2 y 



(B24) 

v 7, J 7 "'' ^ . (B25) 

y/cr q ,+ - °q- V i<l u /v a <i,+ ~ v 1" J \f a q,+ ~ °? - \ -iqu/V v 1" ~ a i - J 

Deviations from the saddle-point loop are expressed in terms of the eigenvectors as 

Sn = (1/VN) ^(c g ,+e,,+ + c 9 ,_e 9 ,_) . (B26) 

Since e* + = e_ g + and e* _ = e_ 9i _, one finds c* + = c_ 9i+ and c* _ = c_ 9 ._. We can choose co.±, as well as Rec 9i ± 
and lmc q} ± with g = 1, . . . , (JV — 1)/2 (for odd AT), as 2 A real normal coordinates. The Jacobian of the transformation 
to the normal coordinates is 2 Ar_1 . 

The operator a has a negative mode (00,- = — 4fco = — 27rv(l — u 2 /v 2 )) related to the longitudinal deviation as 
cq,- = 5Q\f~N /2. Three zero modes 



e ±1 ,_ = (l/V^+^)e ±l ^ T V m ) and e 0)+ = ( J ) 



(B27) 



correspond to solid rotations of the saddle-point loop around the vectors ej by angles 4n, respectively (i = 1,2,3), 
which can be identified with the normal coordinates as 

4>i = - (2N- 1/2 v/Vv 2 + u 2 ) had,- , $2 = (2JV- 1 / 2 t-/V" 2 + « 2 ) Reci,_ , (j> 3 = c q ,+/(Vn sin(6> /2)) . (B28) 



Finally, after excluding the integration over the negative mode by introducing the (5-function J(2co,_/v A) in the full 
normal-mode integral, we obtain the integral over the transverse deviations from the saddle-point configuration in the 
following form: 

/ dfi 9o {g) [ PCcxp(- K - 1 (C,o-C)) = 2TT N N N n N - 2 ir 3 (v 2 + u 2 )yjv 2 - u 2 /VdetZ, (B29) 

JSO(3) J 

where the determinant deter includes 2 (AT — 2) positive eigenvalues of cr, or cquivalently, all eigenvalues of W(9q). 
Now we are in a position to calculate the ratio of the integrals: 

Iso(3) d ^ (g) I V C ex P (-*r 1 (C, <7fl) . 2 

w ? _ — = V dct 7/ det a 4tt k) B30 

J s2 d/j{n)J X>£cxp(-K HCjT?)) 
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with 

00= / <Wff) = 2tt 2 (1 + u 2 /v 2 )y/l - "7« 2 , (B31) 

JSO(3) 

which, in contrast to the individual integrals and determinants, is iV-independent in the N — > oo limit. Indeed, we 
note that the knowledge of the lowest eigenvalues with q <C N, specified in Eqs. (|B17|) and (|B24|) . allows us to 
calculate the ratio of the determinants by re-grouping the factors 

(N-l)/2 

v /det t T/det 7 = (l/<7i,+) [] 7«7«/K+°«.-) ■ (B32) 

The transformation is legitimate because 7 ? 7 g / (<t 9i +o" 9i _) approaches unity sufficiently fast as q increases if q <C N 
and N — * oo. Then the upper limit of the product can be extended to infinity. In the meantime, the Weierstrass 
factorization theorem [l8| for the representation of an entire function as an infinite product leads to the following 
identity 

oo 

JJ(1 - z 2 /q 2 ) = {nz)- l (l - z 2 )- 1 sin^z . (B33) 

q=2 

Therefore, applying this formula with z = u/v and z — ► 1 to Eq. (|B32[) . one finds 

oo 

Vdetcr/ det 7 = (27r) _1 v (v 2 + u 2 )' 1 u 2 v~ 2 q- 2 )/(l - q~ 2 ) = tTVu -1 ^ 4 - w 4 )" 1 sm(Tru/v) . (B34) 

q=2 



APPENDIX C: TWO-CHANNEL SINGLE-STATE MARKOV CHAIN 



In this Appendix we re-derive the statistics of the current working with the simple stochastic model consisting only 
of one site linked to itself by an edge, with jump rates k+ and k_ in the positive and negative directions, respectively. 
This simple model, represented in Fig. [4l describes a continuous time random walk [371 ] . and reflects the instanton 
mechanism of the current generation of our spin-chain model in the weak-noise limit. 

A trajectory in the model is specified by a sequence of the jump directions Wj = ±1 and the times when they occur: 
77 = {u>i, n; w 2l t 2 ; . . . ;w„,t„}, where < n < r 2 < . . . < t. 

To study the statistics of the jump rate one introduces counting stochastic process |33| c(ry) that measures the 
difference of jump counts in the positive and negative directions (which is similar to the definition in Eq. IJT])). Its 
value increases by Wj = ±1 at the jump j if a jump occurs in the positive or negative direction, respectively: 

cfa) = 5>i- (CI) 

j 

We are interested in the statistics of the empirical current tu, i.e., the number of jumps per unit time c(rf)/t. The 
probability of having ut jumps during time t is given by the discrete distribution 

i>,t) - (<W»)> = j ^ e " lnZ < e-c(T,)lnZ > > ( C2 ) 

where S is the Kronecker symbol, and the angular brackets denote the average over stochastic trajectories. This 
average is determined by the Markovian measure Prjexp (— S(r))), where V-q = J2 n ■ • • anc ^ the action reads 

n 

exp(-5( J7 )) = n^ e " (K++ ^ )t ' ( C3 ) 

3=1 

with k±i = k± being the jump rates. 

The genera ting function [e~ c ^ ln z ) of the distribution P(w,t) can be calculated using the standard procedure, e.g. 
described in [311 ]: 

(e- £( "» lnZ ) =e x ^\ (C4) 
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where 

\{Z) = Z~ 1 k+ + Zk_ - k+- k_ (C5) 

is the only eigenvalue of the biased lxl transition "matrix", obtained by replacing the jump rates as n + — ► Z~ 1 k + 
and k_ — * while keeping the overall escape rate equal to {k+ + 

Although the distribution can be obtained exactly in terms of the modified Bessel function, 

P(uj, t) = I ut {2jK^Zt) (K+/^r t/2 e -( K + + K -)* , (C6) 

to derive the Cramer function S(u>) = lim^oo t _1 In P(u>, t) we only need to calculate the integral in Eq. (|C2p within 
the saddle-point approximation. Thus, one derives 

S((j) = -X(Z)-(j]nZ, (C7) 

where Z is expressed in terms of u> by means of the saddle-point equation 

d z S = -dzX(Z) - luZ- 1 = (C8) 

with X(Z) specified in Eq. (|C5[) . The resulting quadratic equation for Z has two real solutions of opposite signs. The 
positive solution should be chosen 



Z = V^ 2 + ^--^ (C9 ) 
2k— 

since it corresponds to the minimum of S over the integration contour. 

The expressions presented above provide a link between the spin-chain model in the weak noise limit and the 
single-state Markov chain. 
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